Non exponential quasiparticle decay and phase relaxation 
in low dimensional conductors 
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We show that in low dimensional disordered conductors, the quasiparticle decay and the relaxation 
of the phase are not exponential processes. In the quasi-one dimensional case, both behave at small 
time as e _< - i/ ' Ti "- ) ' 1 where the inelastic time Ti n , identical for both processes, is a power T 2//3 of 
the temperature. The non exponential quasiparticle decay results from a modified derivation of the 
Fermi golden rule. This result implies the existence of an unusual distribution of relaxation times. 
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The issue of dephasing in the presence of electron- 
electron interactions in disordered conductors is of great 
importance in mesoscopic physics. This problem, first 
addressed by Altshuler, Aronov and Khmelnitskii (AAK) 
[1], has been recently revisited in the light of a new set 
of experiments [2-4]. A related problem is the under- 
standing of the time evolution of a quasiparticle state 
due to electron-electron interaction, which governs the 
relaxation towards thermal equilibrium [5]. Since coher- 
ent effects in disordered systems result from the coherent 
pairing of two scattering amplitudes defined for a given 
quasiparticle state, the coherence is lost once this state 
has relaxed. Thus it seems natural to assume that quasi- 
particle and phase relaxations are of the same nature 
[6] and are driven by the same time scale. In this let- 
ter, we show that in low dimensional disordered conduc- 
tors and in particular for quasi-one-dimensional (quasi- 
ld) wires, both relaxations are faster than exponential 
and are driven by the same characteristic time [7]. This 
non exponential behavior reflects the existence of a dis- 
tribution of relaxation times. Such a non exponential de- 
cay is unusual in quantum condensed matter physics but 
more frequent in the context of molecular relaxation pro- 
cesses and in glassy systems. Stretched and compressed 
exponentials are mostly used as a way to fit unusual re- 
laxations but no microscopic basis can be assigned to 
account for this behavior [8]. Here, we derive it from a 
new treatment of the Fermi golden rule. 

We shall first consider the quasiparticle decay, using 
the Fermi golden rule [9] . We show that, due to screened 
Coulomb interactions with small energy transfer, the re- 
laxation rate is not constant, implying a non exponen- 
tial decay. This results from a key step in the Fermi 
golden rule which stems that the transitions conserve en- 
ergy within h/t where t is the duration of the perturba- 
tion. Usually this constraint is of no practical importance 
and energy conservation is described by a delta function. 
Here we show that this is no longer possible. As a re- 
sult, we find that for quasi-ld wires, the probability for 
a quasiparticle to stay in its initial state behaves, at small 



times, as V(t) = e - ^*/ 7 "*") 3 . The temperature depen- 
dence of the inelastic time is Ti n (T) cx T~ 2 / 3 and (3 is a 
numerical constant. 

Then, we shall come to the relaxation of phase coher- 
ence. It is given by the average (e**^) of the relative 
phase between time reversed trajectories. Starting 
from the AAK calculation, we show that the phase relax- 
ation is also non-exponential and that, at small times, it 
behaves like (e 1 *'*- 1 ) ~ e -( t / T 4-f / where the phase coher- 
ence time Tcf, is proportional to the quasiparticle decay 
time T in . 

We start by considering the decay of a quasiparticle, 
recalling first some known features of the derivation. Us- 
ing the Fermi golden rule, the quasiparticle lifetime can 
be written in terms of a kernel K(u>) which is the aver- 
age over disorder of the squared matrix element of the 
screened Coulomb interaction [7, 10]: 

— - = 8w* cL,K(u,)——- , (1) 
Tin{T) Jo smhphu 

where vq is the total density of states per spin direction. 
The temperature dependence results from the occupancy 
of the initial and final quasiparticle states. Upon disorder 
averaging, the kernel K(cu) is obtained as the squared 
product of the dynamically screened interaction and of a 
long range contribution called diffuson. As a result, we 
have [5] : 

The diffusive nature of the electronic motion implies a 
strong dependence of the transition probability upon the 
space dimensionality d that appears in the sum over 
the modes q. The kernel K(ui) then depends on the 
space dimensionality d and it is given by [5] K(lo) — 
(a d /l6^Lu 2 ) (uj/E c ) d/2 where a x = \/2/tt 2 , a 2 = 1/2tt 2 , 
(X3 = \/2/27r 3 and E c = D/L 2 is the Thouless frequency. 
For d — 3, the integral in (1) is convergent so that Ti n {T) 
is well defined and behaves like T~ 3 / 2 . However for d < 2, 
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the integral in (1) diverges at low energy transfer u>. To 
cure this divergence, it is commonly argued that the low 
frequency cut-off needed is l/r in itself, since no energy 
transfer can be smaller than h/ri n . Consequently, the 
lifetime is solution of a self-consistent equation whose so- 
lution in d = 1 is : 



1 



Sah 2 



2/3 



(3) 



where the conductivity is a = 2e 2 voD / {LS). This tem- 
perature dependence has been first obtained by Altshuler 
and Aronov [5]. 

We argue here that this divergence is indeed the sig- 
nature of a new behavior for the quasiparticle decay. We 
prove that this decay is actually non exponential. The 
crucial point in our argument is that it is not correct to 
cut-off the integral at 1/T in . To grasp the relevance of 
this statement, it is important to recall the Fermi golden 
rule prescription namely that, after a time t, the range 
of accessible states is limited to energies larger than h/t, 
not h/Tm- This leads to the replacement of Eq.(l) by the 
following expression for the disorder averaged transition 
probability V^ 2 \t) towards final states, calculated up to 
second-order in perturbation : 



(9) = tt^t r 

2vtfi 2 h 
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duj 
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d/2 



(4) 



where, for simplicity, the thermal factor has been re- 
placed by an upper cut-off at hui <~ T. We have used 
the above expressions for the kernel in d dimensions. In 
one dimension, this leads immediately to a t 3 ! 2 power 
law : V (2) {t) = — ^=t 3 / 2 so that the quasiparticle 
relaxation is not exponential. 

Let us prove now this qualitative statement, coming 
back to the derivation of the Fermi golden rule. A given 
initial quasiparticle state a interacts with a quasiparticle 
of energy e^, leading to two quasiparticles of final energies 
ep and eg. As known from quantum mechanics textbooks 
[11], the transition probability towards final states is : 



v( 2 Ht) = 2j2\u cn ,ps\ 2 M 

Pit 



e Q + e 7 — £/3 — £,5 



) (5) 



where U aj jjs is the matrix element of the interaction. 
The function / t (Aw) of width w/t is given by [11] 



A (Aw) = 



sinAwt/2\ 2 
Aw/2 ) 



(6) 



Its maximum is equal to t 2 and its integral is 2irt. Usu- 
ally, this function can be approximated by 27rf 5(Au>), 
so that the decay is linear in t, and the prefactor is the 
inverse quasiparticle time. 

The main idea here is that this approximation is not 
always valid. To see this, we first calculate the disorder 




FIG. 1: Plot of the function h(ut). It vanishes for small 
argument, justifying the cut-off of order 1/t in Eq. (4). 



average of Eq. (5) using standard methods [7]. The new 
input is that the energy of the initial and final states 
may differ by a small amount Auo of order 1/t, as ex- 
plicitly seen on Eqs.(5) and (6). As a result, the diffu- 
sons and dynamically screened interactions which enter 
in the kernel K have to be taken at different frequen- 
cies uj± = uj ± Aw/2. This immediately leads to a kernel 
i^Aufw) which now depends on the energy difference Aw : 
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4tt 2 ^ 

q^o 



yj (w 2 + D 2 q i )(uj 2 _ + D 2 q i ) 



(7) 



instead of (2). This kernel yields to a transition proba- 
bility of the form : 



v (2 Ht) 



T/h 



dujg t (uj) 



where 



9t(oj) 
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(8) 



(9) 



Letting ft(Au) — 2irt5(Acu) leads immediately to the 
usual result, namely a behavior of p( 2 \t) linear in t and 
defines l/r in provided the integral is convergent. 

However in Id, when Auo = 0, K AuJ=0 {u) = K(u>) oc 
1/uj 3 / 2 and the integral on u becomes divergent. It is 
thus crucial to keep the full expression of f t (Aw). Doing 
this, we find that for u ^> 1/t, we still have gt(ui) = 
2ntK(u>) oc t/uj 3 / 2 , while for uj <C 1/t, it is easy to check, 
since /*(0) = t 2 , that the function gt{oj) oc t 2 /^/uo, so 
that its integral near zero frequency is indeed convergent. 
More precisely, we can show that gt(ui) is of the form 
<7t(w) = 2-KtK(u))h(ujt) where the function h(u>t) is linear 
for small argument and tends to 1 for large argument. 
This function is calculated numerically and is shown on 
Fig. 1. The quasiparticle decay probability takes the 
form 



p(2)( t ) 



Tt 



irV2h 2 



duo 

3/2 



h(ut) 



(10) 



We have thus proven that the cut-off in (4) appears 
as a natural consequence of a proper use of the Fermi 
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golden rule. More precisely, because of the function h 
which naturally provides a lower cut-off of order 1/t, the 
integral now converges and 

At small times, the survival probability V(t), i.e. the 
probability that a quasiparticle stays in its original state, 
is given by V(t) = 1 — V^ 2 '(t). By exponentiating this 
relation, we obtain : 

V(t) = er f3 (^t) 3/2 . (12) 

The survival probability is thus given by a compressed ex- 
ponential characterized by the inelastic time Eq.(3) and 
where j3 ~ 5.83\/2/7r. Similarly, in two dimensions for 
a film of thickness a, using (4), one finds a logarithmic 
correction V{t) ~ e^^t^ with t^-{T) cx These 
temporal behaviors constitute one of the main results of 
this paper. We emphasize again that there are direct 
consequences of the Fermi golden rule prescription ac- 
cording to which the energy is conserved, not within the 
decay rate fr/T in , but rather within the inverse time h/t 

[11]- 

Now, one can wonder whether this peculiar behavior of 
the energy relaxation has its signature in the time depen- 
dence of the phase relaxation of coherent effects in weakly 
disordered systems. These effects result from the coher- 
ent pairing of two scattering amplitudes defined for a 
given quasiparticle state. In particular, we consider pairs 
of time reversed trajectories (the cooperon) as it appears 
in the weak-localization correction to the conductivity. 
The relaxation of the cooperon is driven by a character- 
istic time called the phase coherence time. It seems 
quite intuitive that, as far as Coulomb interactions are 
involved, the quasiparticle decay and the cooperon relax- 
ation are related. Thus the question arises of the relation 
between the inelastic time r in and the phase coherence 
time T0. We shall now show that these two relaxation 
processes are indeed identical and characterized by the 
same time scale. 

To that purpose, we consider the time relaxation of 
the cooperon by replacing the Coulomb interaction by a 
classical fluctuating potential V(r, t) whose characteris- 
tics are determined by the fluctuation-dissipation theo- 
rem [1]. The cooperon contribution to the return proba- 
bility can be written under the form : 

P c (r,r,t) = Pi°)(r,r,i)(e^) Tc (13) 

where Pi°' is the cooperon in the absence of the fluctu- 
ating potential and <& = $(r, t) is the relative phase of a 
pair of time reversed trajectories at time t : 

^ = ^j\v(r(r),T)-V(r(r),r)}dr (14) 



This expression is valid in the eikonal approximation i.e. 
for a slowly varying potential whose effect is to multiply 
the disorder averaged Green function by a phase term 
proportional to the circulation of V(r(r),r) along the 
trajectory between the times and t. We define r = t — t. 

The dephasing $ is accumulated along the diffusive 
electronic trajectories paired in the cooperon. One of 
them propagates in the time interval < t < t whereas 
its time reversed counterpart propagates from r = t to 
r = 0. We denote by (• • -}t,c the average taken both over 
the distribution of the diffusive trajectories ((• • -)c) and 
over the thermal fluctuations ((• • -)t) of the electric po- 
tential. The latter are Gaussian so that the thermal aver- 
age (e J *) T = e~^(* ) T . Using (14) and the fluctuation- 
dissipation theorem in the classical limit (f3huj <C 1), 
namely, (l^(q, cj)V(— q, w)) T = , we obtain for the 
dephasing the following expression 

4g2T C ^ f dq l-cosq.(r(T)-r(r)) 

(15) 

The average (cosq.(r(r) — r(r))) over the diffusive tra- 
jectories of time t is e~ 2Dq T \ 1 - 2r / t \, F or a quasi-ld wire, 
the integrations over q and r lead to 

Assuming first that ^e~ 5 (* ~ e 2 ^* ^ T - c , we ob- 

tain for the cooperon, at small time, the compressed ex- 
ponential behavior 

( e **) TC ~ e -V^/4(t/n„) 3 / 2 (17) 

identical to the energy relaxation (12) and with the 
same characteristic time Ti„ given by (3). A similar be- 
havior for the phase relaxation has been also found in 
[12]. It is interesting at this stage to compare (15) with 
(4) obtained for the transition probability of a quasi- 
particle state. Although these expressions behave sim- 
ilarly, the convergence in Eq.(4) results from a cut-off 
at small ui whose origin is in the Fermi golden rule pre- 
scription namely, that among the large number of ac- 
cessible states in d = 1, only those with energy trans- 
fer larger than h/t are accessible after a time t. This 
low energy cutoff does not exist in (15) and the con- 
vergence results from the compensation between the two 
terms in the bracket that describe respectively the contri- 
butions of the correlations (V(t(t),t)V(t(t'),t'))t and 
(V(r(r),r)V(r(r'),r / ))T to the cooperon. 

The result (17) is not fully correct since we have ap- 
proximated the average (• • -)c of the exponential by the 
exponential of the average. Using the functional inte- 
gral approach presented in [1], there is a way to de- 
rive an expression for the phase relaxation valid at all 
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FIG. 2: Behaviour of (e l<s< -^) T ,c ■ The continuous line is the 
exact result (20). The dotted line is obtained from the small 
time expansion (17). The dashed line shows the exponential 



times by considering the Laplace transform P 7 (r, r) = 
J dtP c (r, r, t)e~ 7 * of the cooperon. In quasi-ld, one has : 



P 7 (r,r) - 



Ai (r in /r 7 ) 

D Ai'(T in /T 7 ) 



(18) 



where Ai et Ai' are respectively the Airy function and its 
derivative [13] and t 7 = I/7. The probability P c (r,r,t) 
in (13) can thus be obtained from the inverse Laplace 
transform of (18). Since, in the quasi one-dimensional 



limit, one has Pi \r,r,t) = l/v^wDt, the dcphasing 
term (e l *) TC is a function of t/r,„ that satisfies 



dt 



'T,C 



-t/T~, _ _ 



= -\fTtTi; 



. Ai(r^ w /r 7 ) 
Ai' (r in /r 7 ) 



(19) 



The inverse Laplace transform is obtained by noticing 
that the Airy function and its derivative are analytic and 
non meromorphic functions whose zeroes lie on the neg- 
ative real axis. Then, by performing the integral in the 
complex plane with the residues Res(e s *Ai(s)/Ai'(s)) = 
e~l" n l*/|u„| where the u n are the zeros of Ai'(s) given 

at a very good approximation by \u n \ — (^-(n— §)) 2 ^ 3 
[13], we obtain the analytic function 



-kr.lt/Ti, 



(20) 



n=l 



At small times t < Tj n , it behaves like Eq.(17). At large 
time, the relaxation is driven by the first zero of the 
Ai' function, namely <e l *) TC ~ ^/nt/T^e-^ 1 /^ /\m\ 
with |ui| ~ 1.019. Clearly, the relaxation (20) is never 
exponential. It appears as a distribution of relaxation 
times Ti n /\u n \ which is at the origin of the rather un- 
expected compressed exponential behavior of the quasi- 
particle decay and of the cooperon phase relaxation. The 
expression (20) constitutes one of the main results of this 
paper. The question arises of how this behavior could 
show up experimentally. It has been stressed by Pierre 



et al. [4], that the Laplace transform (18) is well approxi- 

r \ 1/2 

mated by the relation P 7 (r, r) = (1/2-/D) ( ^ + 

where a ~ 1/2 is an adjustable numerical constant. This 
approximation corresponds to an exponential relaxation 



-t/2r«, 



that is clearly at odd with the be- 



havior (20) (see Fig. 2). However, the difference be- 
tween the exact relation and the exponential approxi- 
mation is difficult to sec experimentally The time r 7 
accounts for other processes such as the decay rate in 
a magnetic field which, for a wire of section S, is given 
by r" 1 = t b 1 = DS 2 e 2 B 2 /3h 2 [14]. A possibility to 
probe the i 3 / 2 behavior at small time is to study the 
limit tb <C T in where the weak localization correction to 
the conductivity is 



At = — \/~Dtb 

nn 



1 (TB 

4 \ Ti„ 



3/2' 



(21) 



The power-law dependence of the correction term in (21) 
is a direct signature of the t z l 2 behavior of the relax- 
ation at small time. The asymptotic behavior Acr(P) — 
Acr(P -> 00) cx T/B 4 reflects both the T~ 2 / 3 dependence 
of Tin and the non exponential i 3 / 2 phase relaxation. An- 
other direct probe of the non exponential relaxation of 
the phase should may also be provided by the behavior 
of the ac conductivity a(u). 
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